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Phylogenetic Analysis of Feline Coronavirus Strains in an Epizootic 
Outbreak of Feline Infectious Peritonitis 

E.N. Barker, S. Tasker, T.J. Gruffydd-Jones, C.K. Tuplin, K. Burton, E. Porter, M.J. Day, R. Harley, 

D. Fews, C.R. Helps, and S.G. Siddell 

Background: Feline coronavirus (FCoV) infection is common. In a small percentage of cats, FCoV infection is 
associated with the fatal disease feline infectious peritonitis (FIP). Genetically distinct virulent and avirulent strains of 
FCoV might coexist within a cat population. 

Objectives: To determine whether the strains of FCoV in FlP-aflfected cats are closely related or genetically distinct 
from the fecally derived strains of FCoV in contemporary-asymptomatic cats during an epizootic outbreak of FIP. 

Animals: Four cats euthanized because of FIP and 16 asymptomatic cats. 

Methods: This prospective outbreak investigation was initiated during an outbreak of FIP in cats within or rehomed 
from a rescue/rehoming center. Postmortem samples were collected from cats with FIP and contemporaneous fecal samples 
from asymptomatic cats. RNA was purified from tissue and fecal samples, FCoV gene fragments were reverse transcribed, 
PCR-amplified using novel primers, and sequenced. Sequences were aligned with ClustalW and compared with published 
FCoV sequences. 

Results: FCoV RNA was detected in all 4 FIP cat postmortem samples and in 9 of the 16 fecal samples from con- 
temporary-asymptomatic cats. Novel primers successfully amplified fragments from 4 regions of the genome for all 
FCoV-positive samples. Phylogenetic analysis showed that the FlP-associated strains of FCoV from the outbreak were 
very closely related to the fecally derived strains of FCoV from contemporary-asymptomatic cats. 

Conclusions and Clinical Importance: Sequence analysis provided no evidence that genetically distinct virulent and avir¬ 
ulent strains of FCoV were present during this FIP outbreak. 

Key words: Infectious disease: Polymerase chain reaction; Sequence analyses. 


F eline coronavirus (FCoV) infection is endemic 
among cats worldwide. In the United Kingdom, 
around 40% of the domestic cat population has been 
infected with FCoV and in multicat households, this 
figure increases to 90%. 1-3 Natural infections with 
FCoV are often transient and asymptomatic or result 
in mild gastrointestinal disease. In these cases, the 
causative FCoV is sometimes referred to as feline 
enteric coronavirus. However, in some animals, FCoV 
infection results in the systemic disease of feline infec¬ 
tious peritonitis (FIP), 4 and in these cases the virus 
has been referred to as FIP virus. FIP usually arises 
sporadically and unpredictably within a subpopulation 
or a group of cats. 4 Epizootic FIP, in which FIP arises 
at an above average frequency in a single locality over 
a short period of time, is far less common. 4 ' 5 

The pathogenesis of FIP has been investigated for 
over 40 years, 4 but remains incompletely understood. 
Viral factors are thought to be important: in experi¬ 
mental studies FCoV obtained from FlP-affected cats 
and FCoV obtained from the feces of asymptomatic 
cats manifest their respective pathogenic phenotypes 
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(FIP or no to mild enteric disease, respectively) when 
inoculated into cats. 4 Accordingly, 1 recent hypothesis 
states that there are distinct populations of virulent 
(FlP-associated) and avirulent (largely localized to the 
intestines) FCoV strains circulating in cat populations, 
and that these are independently acquired: this is 
known as “the circulating virulent/avirulent” hypothe¬ 
sis. 6 An alternative hypothesis states that FIP arises 
within individual infected cats as a result of mutation 
(s) in the virus genome: this is known as “the internal 
mutation” hypothesis. 7-10 However, no specific viral 
mutation has been identified that is associated with the 
FIP phenotype. 

Phylogenetic studies supporting each of these hypoth¬ 
eses have been published. 6 ' 8 ’ 10 The circulating virulent/ 
avirulent hypothesis was supported by the phylogenetic 
clustering of membrane glycoprotein and 7b gene frag¬ 
ment sequences from tissue-derived FCoV strains from 
FlP-affected cats (n = 8) and fecally derived FCoV 
strains from asymptomatic cats (n = 48). 6 In contrast, 
other studies analyzing sequence data from membrane 
glycoprotein 8 and accessory protein 3c 11 genome frag¬ 
ments from FCoV strains collected from the feces of 
asymptomatic cats (n = 15 and 27, respectively), and 
from the tissues of FlP-affected cats (n = 28 and 28, 
respectively), found phylogenetic clustering according to 
the geographic location of sampling rather than in asso¬ 
ciation with the disease phenotype. However, both of 
these studies are difficult to interpret because sampling 
occurred over protracted periods of time (2-3 years and 
unspecified), at different geographic locations, and from 
cats with different genetic and environmental back¬ 
grounds. 

Recently, the Feline Centre, Langford Veterinary 
Services, University of Bristol was able to obtain 
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samples from an epizootic outbreak of FIP at a single 
site cat rescue and rehoming center in the UK. Sam¬ 
pling involved kittens entering the center from July 
2010 to January 2011. These samples provided the 
unique opportunity to investigate the phylogeny of the 
FCoVs present in an epizootic outbreak of FIP, which, 
to our knowledge, has not been reported previously. 
Thus, the aim of this study was to determine whether 
the strains of FCoV in FlP-affected cats were closely 
related or genetically distinct from the fecally derived 
strains of FCoV in contemporary-asymptomatic cats 
during this FIP outbreak. 

Materials and Methods 

Samples and RNA Purification 

Samples were obtained from 4 kittens with suspected FIP. 
Postmortem examination was performed immediately after eutha¬ 
nasia of 3 kittens with FIP (named Z/FIP, F/FIP, and J/FIP; 
euthanized January 19, January 29, and February 4, 2011, respec¬ 
tively), at which time tissue samples were collected and stored in 
10% neutral buffered formalin for routine histopathology, and 
RNA later 3 (according to manufacturer’s instructions) for RNA 
purification. Ascitic fluid collected as part of a diagnostic investi¬ 
gation was available for another kitten with suspected FIP 
(U/FIP; euthanized January 30, 2011). In January 2011, fecal 
samples were collected from 16 clinically healthy kittens and cats 
(identified by number) at the same rescue and rehoming center 
and stored immediately at —80°C. Total RNA was purified from 
20 mg of tissue (omentum from kitten J/FIP; colonic lymph node 
from kitten Z/FIP; mesentery from kitten F/FIP), 100 pL of as¬ 
citic fluid (from kitten U/FIP) and from 0.5 g feces (clinically 
asymptomatic kittens) using a commercially available kit (Nucleo- 
Spin RNA II b ) as described previously. 12 Formalin fixed tissue 
from the 3 kittens with FIP for which samples were available 
(Z/FIP, F/FIP and J/FIP) were subjected to immunohistologic 
examination for the presence of FCoV antigen. 13 All kittens had 
been resident at a single site cat rescue and rehoming center in 
the UK. The following pairs of kittens were litter mates and resi¬ 
dent within the same pen: kittens Z-FIP and #80; #65 and #73; 
#66 and #75. Otherwise the kittens were housed in separate pens 
within the same complex. Each litter of kittens had been relin¬ 
quished to the center from different private households and kit¬ 


tens that subsequently went on to develop FIP had been 
rehomed separately to private households before the onset of 
clinical signs. 

Reverse Transcription PCR and Sequencing 

FCoV RNA was amplified by real-time, reverse transcriptase 
polymerase chain reaction (RT-qPCR) as described previously. 12 

To design gene specific (GSP) reverse transcription and PCR 
primers, all of the FCoV genome sequences found on the 
SARGENS website http://veb.lumc.nl/SARGENS/ accessed May 
28, 2011 were aligned using ClustalW in MacVector vl2° and 
primers (Table 1) were designed to 4 target regions within the 
FCoV genome (Fig 1) using Primer 3. 14 The targets comprised of 
the following: 2 regions within the replicase polyprotein; an 
amino-proximal region of nsp3, and the carboxyl-proximal region 
of nsp 12 (RNA-dependent RNA polymerase); an amino-proximal 
region of the surface glycoprotein; and the membrane glycopro¬ 
tein. All primers were synthesized by Metabion International 
AG d . 

Purified RNA from each sample was reverse transcribed 
using the ImProm-II Reverse Transcription System 6 using 
250 ng random hexamers or 500 nM GSP reverse primer, 
500 pM each dNTP. and 1-9 pL RNA combined with water to 
a total volume of 20 pL according to manufacturer's instruc¬ 
tions. The resultant cDNA was amplified by PCR using the 
iProof High-Fidelity PCR Kit 1 using 500 nM of each primer, 
500 pM each dNTP. and 1 pL cDNA template, combined with 
water to a total volume of 50 pL. The following thermal cycling 
protocol was used for the nsp3, nsp 12, and the membrane 
glycoprotein gene fragments: 98°C for 30 seconds and 40 cycles 
of 98°C for 10 seconds, 59°C for 20 seconds, and 72°C for 
30 seconds. The following thermal cycling protocol was used for 
the surface protein gene fragment: 98°C for 30 seconds and 40 
cycles of 98°C for 10 seconds, 62°C for 20 seconds, and 72°C 
for 45 seconds. After a failure to amplify the surface protein 
gene fragment from the Z/FIP sample using iProof polymerase, 
the HotStarTaq Master Mix Kit 8 was used with 200 nM of pri¬ 
mer (FCoV SI Fl/Rl), and 1 pL cDNA template in a total 
volume of 25 pL, with cycling conditions of 95°C for 15 min¬ 
utes and 45 cycles of 95°C for 10 seconds, 56°C for 20 seconds, 
and 72°C for 120 seconds. Positive (FCoV cDNA) and negative 
(H 2 0) PCR controls were included in each reaction. Reactions 
were carried out in either a PTC 200 DNA Engine 11 or MJ Mini 
thermal cycler 1 . 


Table 1. Primers used in directed reverse transcription and PCR reactions. 


Target Site on Genome 

Direction 

Primer Name 

Oligonucleotide Sequence 

Estimated 
Product Size 3 

Replicase polyprotein (amino 

Forward 

FCoV nsp3 A FI 

5'-ATCCATATGGTTCTGGCATGG-3' 

730-970 bp 

-proximal region of nsp3) 

Reverse 

FCoV nsp3 A R2 

5'-TTTAGCYGTACTATAATCATTGAGCA-3' 


Replicase polyprotein (carboxyl 

Forward 

FCoV nsp 12 B FI 

5'-CCCACAATGACTCAAATGAA-3' 

800 bp 

-proximal region of nsp 12) 

Reverse 

FCoV nsp 12 B R1 

5'-TCTGGTTCYACCCAACACTT-3' 


Amino-proximal region of 

Forward 

FCoV SI FI 

5'-TCTGTKGCCATCAAAATCAC-3' 

1900 bp 

the surface glycoprotein 11 

Reverse 

FCoV SI R1 

5'-CATTAACATCHACCATTACATCTG-3' 



Forward 

FCoV SI FB 

5'-GGAAGAGAATCAGCCTCACG-3' 

Sequencing 


Forward 

Reverse 

Reverse 

FCoV SI FC 

FCoV SI RB 

FCoV SI RC 

5'-TTGCGCTGGTTATGCTAAGA-3' 
5'-CACGACCCTGT ACCA AT GTG-3 1 
5'-CACCTGTCCCACAGTATGGT-3' 

primers 

Membrane glycoprotein 11 

Forward 

Reverse 

FCoV M FI 

FCoV M R1 

5'-GCGGTTMTAAACGAAATTGA-3' 

5'-TGAGTAATCACCRGCTTTAGATTT-3' 

1040 bp 


“Based on available feline FCoV genome sequences. 

b Sequence variability in the 5' region of the surface and membrane glycoproteins necessitated the placement of the forward primer in 
the preceding highly conserved regions of the nspl6 and small envelope protein genes respectively. 
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Fig 1 . Feline coronavirus genome with component genes, nucleotide scale, and for each of the 4 primer pairs their approximate binding 
sites, their target gene fragment, and the predicted amplicon size. 


PCR amplicons were purified (NucleoSpin Extract II b ), quan¬ 
tified 1 , and sequenced 1 . Because of the length of surface protein 
gene amplicon, additional sequencing primers were designed 
using Primer 3 (Table 1). 

Data Analysis 

Nucleotide sequences were aligned using ClustalW, 15 in Mac- 
Vector vl2, with subsequent manual adjustment. Phylogenetic 
trees were constructed using the neighbor-joining program from 
a distance matrix, 16 corrected for nucleotide substitutions by the 
Kimura two-parameter model. 17 Further phylogenetic trees were 
constructed using maximum likelihood (PhyML v3.0) and parsi¬ 
mony (dnapars v3.67) methods accessed via the Mobyle portal at 
http://mobyle.pasteur.fr/. The data set was resampled 1,000 times 
to generate bootstrap percentages. Sequence data corresponding 
to the same regions of the canine coronavirus genome 
(GQ477367) were used to root the phylogenetic trees. The FCoV 
fragments were translated in silico into their corresponding 
amino acid sequences by Transeq accessed via EMBL at http:// 
www.ebi.ac.uk/Tools/emboss/transeq/, using the standard transla¬ 
tion table, then aligned using ClustalW. The FCoV gene frag¬ 
ment sequences were deposited in the European Molecular 
Biology Laboratory Nucleotide Database (accession numbers 
HE860440-91). 

Ethics 

All of these samples and details were obtained with the full 
consent of the rescue and rehoming center or owners as appropri¬ 
ate. In the FIP cases, euthanasia was performed on the basis of 
veterinary advice. 

Results 

All 3 kittens suspected of having FIP and for which 
tissues were available had histopathologic changes con¬ 
sistent with FIP (pyogranulomatous parenchymal foci, 
perivascular mononuclear infiltrates, fibrinous polyse¬ 
rositis), and immunohistologic evidence of FCoV anti¬ 
gen within tissue-associated macrophages, confirming a 
diagnosis of FIP. FCoV RNA was detected by RT- 
qPCR in the tissues (threshold cycle values (Cts): J/ 
FIP 15.5, F/FIP 17.1, Z/FIP 17.3) from the kittens 
with confirmed FIP and the ascitic fluid (U/FIP Ct 
29.1) from a kitten with suspected FIP. Nine of the 16 
fecal samples from asymptomatic kittens and cats 
(56%; identified as #65, #66, #67, #71, #73, #75, #76, 
#77, and #80) were positive for FCoV on RT-qPCR 


(Ct 18.6-27.8). Of these, four (#66, #71, #76, and #77) 
have been euthanized (see Discussion for details). The 
remainder was recorded as alive on or after March 30, 
2012 . 

Using the primer pair nsp3 A F1/R2, a 896 bp frag¬ 
ment of the nsp3 gene was amplified by PCR for each 
of the FCoV-positive samples, from which sequence 
data were obtained. Nucleotide identities between the 
sequences ranged from 99.6 to 100%. All of the amp¬ 
licon sequences clustered together in a distinct clade 
when aligned with other FCoV genomes (Fig 2A). 

Using the primer pair nspl2 B Fl/Rl, a 803 bp 
fragment of the nspl2 gene was amplified by PCR for 
each of the FCoV-positive samples, from which 
sequence data were obtained. Nucleotide identities 
between the sequences were identical for all but one of 
the samples, which had a single synonymous substitu¬ 
tion at one site. All of the amplicon sequences clus¬ 
tered together in a distinct clade when aligned with 
other FCoV genomes (Fig 2B). 

Using the primer pair SI Fl/Rl, a 1,881 bp frag¬ 
ment comprising a short 3' terminal region of the repli¬ 
case polyprotein (subunit nspl6) gene and the 5' end 
of the surface glycoprotein gene (providing approxi¬ 
mately 32% coverage of the surface protein gene of 
FCoV ClJe [DQ848678]) was amplified by PCR for 
each of the FCoV-positive samples, from which 
sequence data were obtained. Nucleotide identities 
between the sequences ranged from 99.0 to 100%. All 
of the amplicon sequences clustered together in a 
distinct clade when aligned with other FCoV genomes 
(Fig 2C). 

Using the primer pair M Fl/Rl, a 1,043 bp frag¬ 
ment comprising a short 3' terminal region of the small 
envelope protein gene, and the majority the membrane 
glycoprotein gene (approximately 90% coverage of the 
membrane glycoprotein gene of FCoV ClJe) was gen¬ 
erated by PCR for each of the FCoV-positive samples, 
from which sequence data were obtained. Nucleotide 
identities between the sequences ranged from 99.2 to 
100%. All of the amplicon sequences clustered 
together in a distinct clade when aligned with other 
FCoV genomes (Fig 2D). 

For all genome fragments, maximum likelihood and 
parsimony methods produced similar phylogenetic 
trees (with equivalent bootstrap values) to those found 
by the neighbor-joining method (data not shown). 
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When sequence data for the outbreak samples were 
compared with each other there were no nucleotide, or 
predicted amino acid, motifs that could differentiate 
between the fecally derived FCoV strains from asympto¬ 
matic cats and the tissue- or ascitic fluid-derived FCoV 
strains from cats with FIP. 

Discussion 

In the epizootic outbreak described in this study, 
postmortem derived tissues or ascitic fluid were avail¬ 
able from 4 cats of which 3 were confirmed to have 
FIP by histopathology and the remaining cat was 
strongly suspected of having FIP based on clinical fea¬ 
tures. In each case, the diagnosis of FIP is supported 
by the amplification of FCoV RNA from their respec¬ 
tive tissue or ascitic fluid samples. Unfortunately, the 
true incidence of FIP was unknown because of the 
high throughput of cats within the center and minimal 
available follow-up data. Of the 16 fecal samples col¬ 
lected from contemporaneous asymptomatic cats, 9 
(56%) contained FCoV RNA. This is similar to what 
was found in a study that identified FCoV shedding in 
60% of cats 1 week after arriving at a rescue shelter. 18 
More than 1 year postsampling, 5 of the 9 contempo¬ 
raneous asymptomatic cats remain alive and without 
any signs that could be suggestive of FIP over 1 year 
postsampling. We therefore suggest that these isolates 
can, with some certainty, be considered as feline 
enteric coronavirus infections based on their biological 
behavior. However, as noted, four of the contempora¬ 
neous asymptomatic cats were euthanized subse¬ 
quently. Two of these cats were euthanized at the 
rehoming center with nonspecific clinical signs; the 
other 2 cats had developed signs potentially consistent 
with FIP (eg, malaise, abdominal fluid) after rehoming 
to new owners and had been euthanized (postmortem 
samples were not available). Consequently, we are 
unable to assign the samples from these animals to the 
FIP or nonFIP groups. This illustrates the caution 
needed during sample collection but does not invali¬ 
date the conclusions of this study. 

The sporadic and relatively infrequent incidence of 
FIP has always confounded the problem of collecting 
multiple contemporaneous samples from cats with FIP 
and asymptomatic cats living in the same location. In 1 
study of 73 FCoV seropositive multicat households 
FIP was encountered at an incidence rate of 4.8% 
within the first 3 years of entering the household. 19 In 
contrast, an epizootic FIP outbreak is defined as occur¬ 
ring where there is an incidence of FIP greater than 


10% in cats within an establishment, 20 and the higher 
incidence of FIP will clearly facilitate the collection of 
relevant samples. However, very few epizootic FIP out¬ 
breaks have been described in the literature. 5 

The primers designed for this study targeted a 
number of regions within the FCoV genome. The sites 
chosen are representative of regions of low (amino- 
proximal region of the surface protein), intermediate 
(membrane protein), and high (carboxyl-proximal 
region of nspl2) conservation, as well as a region that is 
variable in length (amino-proximal region of nsp3). 
Consequently, taken together, these amplicons provide 
an overall picture of genetic relatedness and are appro¬ 
priate for phylogenetic analysis. The primers in combi¬ 
nation with a high-fidelity DNA polymerase (iProof 
High-Fidelity) were successful in generating amplicons 
of the correct size for all samples, except for the tissue 
of 1 FlP-affected cat (Z-FIP). In this case, using the sur¬ 
face protein gene primer pair, a diffuse band of a lower 
molecular weight PCR product was generated. In this 
sample, an alternative DNA polymerase (HotStarTag) 
was subsequently used successfully to generate an 
amplicon of an appropriate size, sequence data from 
repeated reactions were identical. 

All of the FCoV sequences derived in this study, 
both from tissues with typical FIP lesions and from the 
asymptomatic cat feces, clustered together, and sepa¬ 
rately from all other FCoVs that were analyzed. This 
was true for all of the phylogenetic trees generated for 
each of the genome fragments, (Fig 2A-D). Thus, this 
study provides no evidence for the presence of 2 dis¬ 
tinct FIP tissue-associated and fecally derived strains of 
FCoV in the cats associated with this FIP outbreak. 

This conclusion is in contrast to a previous report 
where membrane glycoprotein and accessory protein 
7b gene fragments from FCoV RNA derived from the 
feces of healthy cats formed clusters distinct from 
sequences derived from FlP-associated FCoV RNA. 6 
However, these results are in agreement with another 
study, which found that accessory protein 3c and 
membrane glycoprotein FCoV gene sequences from 
fecally-derived RNA from healthy cats and from RNA 
derived from tissues with typical FIP lesions from spo¬ 
radic cases of FIP were generally clustered according 
to the cattery from where they originated. 811 Our data 
also lead to the conclusion that the cats in this study 
were likely infected with FCoV from the same source 
upon entering the rescue center; although it cannot be 
ruled out that the cats were exposed to very closely 
related strains of FCoV circulating in the geographic 
region before entering the center. 18 


Fig 2. Phylogenetic analysis of nsp3 (A), nspl2 (B), nspl6 and surface glycoprotein (C), and small envelope protein and membrane gly¬ 
coprotein (D) gene fragments for feline coronavirus strains analyzed in this study (in bold) and published feline coronavirus genome 
sequences (as NC_002306 is identical to AY994055 it was excluded from the analyzes) http://veb.lumc.nl/SARGENS/ accessed 28 May 
2011. Phylogenetic trees were constructed by the neighbor-joining method. Evolutionary distances are to the scales shown (number of 
substitutions per nucleotide). The data set was resampled 1,000 times to generate bootstrap percentage values, and values greater than 
80% are given at Ihe nodes of the tree. GenBank accession numbers are shown for all sequences. Fecal sample derived FCoV sequences 
from cats that were subsequently euthanized are marked with an asterix. CCoV, canine coronavirus; FCoV, Feline coronavirus: CLN, 
colonic lymph node. 
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In addition to proposing the “circulating virulent/ 
avirulent hypothesis”, Brown et al 6 suggested that a 
specific amino acid motif within the membrane glyco¬ 
protein could be used to differentiate between the 
different phenotypes of FCoV. 6 This was challenged 
by a later study that found no association between the 
proposed motif and virulent and avirulent pheno¬ 
types. 8 The data generated in this study show that the 
amino acid motif YVIAL (positions 108, 120, 138, 
163, and 199 based on reference sequence for TGEV 
GenBank no. NP058427) was predicted for all of the 
viruses that we analyzed. This motif arrangement was 
not one of the amino acid combinations (YIIAL asso¬ 
ciated with FIP; HIIVI/HIIVL/HVIAL/YVVAL asso¬ 
ciated with asymptomatic cats; and YIVAL seen in 
association with both FIP and asymptomatic cats), 
suggested in the original article. 6 Also, the YVIAL 
motif has been reported subsequently in both FIP- 
associated FCoVs and fecally-derived FCoVs. 8 More¬ 
over, in this study there were no RNA or amino acid 
motifs identified, in any of the gene fragments, which 
were associated with any specific disease phenotype. It 
has to be noted that this study, as with previous phy¬ 
logenetic analyzes of FCoV genomes, was limited by 
only partial genome coverage (approximately 16%). 
Future study involving whole genome sequencing will 
be needed to identify any specific virus determinant 
associated with the development of FIP. 


Footnotes 

a Ambion Ltd, Huntingdon. UK 
b Macherey-Nagel, ABgene, Epson, UK 
c MacVector Inc, Cambridge, UK 
d Martinsried, Germany 
e Promega UK, Southampton, UK 
1 Bio-Rad Laboratories Ltd. Hemel Hempstead. UK 
8 Qiagen, Crawley, UK 
h MJ Research, Waltham, MA 

1 Quant-iT dsDNA Assay Kit, Invitrogen, Paisley, UK 
J DNA Sequencing & Services, University of Dundee, Dundee, 
UK www.dnaseq.co.uk 
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